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Static properties of a bose-fermi mixture of trapped potassium atoms are studied in terms of 
coupled Gross-Pitaevskii and Thomas-Fermi equations for both repulsive and attractive bose-fermi 
interatomic potentials. Qualitative estimates are given for solutions of the coupled equations, and 
the parameter regions are obtained analytically for the boson-density profile change and for the 
boson/fermion phase separation. Especially, the parameter ratio Rint is found that discriminates 
the region of the large boson-profile change. These estimates are applied for numerical results for 
the potassium atoms and checked their consistency. It is suggested that a small fraction of fermions 
could be trapped without an external potential for the system with an attractive boson-fermion 
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I. INTRODUCTION 



Recent development in the cooling and trapping technique of atoms and the following success in achieving the 
Bose-Einstein condensation of alkali atoms have now opened up a new era in the progress in quantum physics. 
While the condensed system serves as a testing ground for a research in fundamental problems of quantum mechanics, 
it also offers a new example of finite quantum many-body systems such as hadrons, nuclei and microclusters. An 
important characteristic of the trapped alkali atoms is that it is a dilute system of weakly interacting particles and 
' an ideal place to test genuine properties of condensed systems predicted by theories. 

Along with a further progress in the study of Bose-Einstein condensed systems, a similar technique is being extended 
to create a degenerate gas of fermionic atoms, where a number of theoretical studies have been made (J] . To realize such 
a degenerate fermionic system that requires a still lower temperature than bose systems, the technique of sympathetic 
cooling has been investigated Q : the cooling mechanism through the collisions with coexisting cold bose particles in 
a polarized boson-fcrmion mixture, where the fermion-fermion interaction is less effective. 

The mixture of bose and fermi particles is itself an interesting system for a study: the hydrogen-deuterium system 
has been studied already at the early stage of these investigations || . Recently, theoretical studies have been made 
for the trapped boson-fermion system of alkali atoms and proposed many interesting properties: the exotic density 
configurations through the repulsive or attractive boson-fermion interactions [fljq], the decaying processes after the 
J-^ | removal of the confining trap ||, and the metastability of the 7 Li- 6 Li mixture lpj. 

In this paper, we study static properties of a mixed bose-fermi system of trapped potassium atoms, where one 
fermionic ( 40 K) and two bosonic isotopes ( 39 K and 41 K) are known as candidates for a realization of such a degenerate 
quantum system. Trapping of the fermionic isotope 40 K has already been reported |0]. These isotopes are also of an 



interest since the boson-fermion interaction may be repulsive or attractive depending on their choice 12 13[ . 

Below we first describe a set of equations for the mixed system at T = 0, where the bosonic part is given by a 
Gross-Pitaevskii (GP) equation and including boson-fermion interaction, and the fermionic part is by the Thomas- 
Fermi (TF) equation. The use of the GP equation allows one to study a system with a net attractive potential for 
bosons where the Thomas-Fermi-type (TF) approximation is not applicable. The applicability of these equations for 
the trapped potassium system was discussed in M . 

In sec. Ill, we discuss qualitative properties of solutions from the GP and TF equations analytically; especially, 
their dependence on the boson/fermion number is estimated. 

In ref. ||, general features of solution including a deformed system is discussed in the repulsive boson-fermion 
interaction. On the other hand, to study more realistic system, we concentrate on the specific systems of potassium 
isotopes, i.e., 39 K- 40 K system and 41 K- 40 K system using empirically derived atomic interactions. We also consider the 
effect of changing parameters of the trapped potential. Finally, we mention the possibility of trapping fermions only 
through the attractive interactions by the trapped bose particles. 



II. FORMULATION 



We consider a spin-polarized system of bosons and fermions at T = described by a Hamiltonian H — Ht+H 
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where </>(r) and V'(^) respectively denote boson and fermion field operators with masses m b and m/. Because of 
diluteness, the particle interactions have been approximated by the s- wave-dominated contact potential v(f — r') = 
{g or h}5 3 (f — r '). The strengths of the coupling constants g (boson-boson) and h (boson-fermion) are given by 

ink 2 , 2irh 2 

9= a bb , h= a b f, (2) 

m b m b f 

where a bb and a b f are the s-wave scattering lengths for the boson-boson and boson-fermion scatterings and m b f = 
m b mf/(m b + rrif) is a reduced mass. The fermion-fermion interaction has been neglected because of diluteness 
and particle polarization. The chemical potentials \x b and [if in ([la]) and (lb) are determined from the condensed 
boson/fermion numbers N b and Nf through the ground state expectation values: 

N b = (Jdr^(r)<j)(r)}, N f = {Jdr^ (r)V(*0)- ( 3 ) 

In the mean-field approximation at T = 0, the bosons occupy the lowest single-particle state <p b (r), and the energy 
of the system is just given by a replacement of field operator <f)(r) in the Hamiltonian with its expectation value: 

$(f) ee (0(r)) = y/Ntipb(f) (4) 

that describes the order parameter of the Bose-Einstein condensate. (We here neglect quantities of the order 0(N b ' 1 ).) 
In this approximation, the boson density n b (r) is given by |<I>(r)| 2 . 

In the same approximation, the fermion wave function is given by a Slater determinant; the single-particle states are 
determined, e.g., by the Hartree-Fock self-consistent equation. In the present approximation of neglecting fermion- 
fermion interactions, we only have to solve for the single particle states under the effective potential for fermions: 
the trapping potential and the boson-fermion interaction. In the actual calculation, we however adopt a semiclassical 
(TF) description for fermion density, which is known to provide a good approximation as far as the number of atoms 
is sufficiently large Q . 

We thus obtain a set of equations for $ and n/ || : 
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where the boson and fermion densities are defined by n b (r) = (0 T (r)</>(r)) ~ |$(r)| 2 and nt(j) = '(r)^(r)) 



III. QUALITATIVE STUDY OF SOLUTIONS 

Before the study of numerical results for (jjj) and (^|) , we discuss about their solutions qualitatively. 

i) Radius of the boson/fermion distribution 

If we neglect the boson-fermion interaction, we would obtain a much broader distribution for fermions as compared 
with that for bosons because of the exclusion principle. To see this, we consider a free boson/fermion system trapped 
in a harmonic oscillator potential. The mean square radius of each system is obtained from the virial theorem: 

(r 2 ) b /N b = \e b , (r 2 ) f /N f = l{n F + 2$, (7) 
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where the harmonic-oscillator lengths £ b j are 



&,/ s W • (8) 

and N b j are the total boson/fermion numbers. The np in (^) is the number of fermions in the last occupied shell that 
is obtained as a solution of Nf = (np + l)(np + 2)(np + 3)/6 ~ np. The ratio of the root-mean-square (rms) radii 

\J{f 2 )f/ { r2 )b t nus increases as the power N 1 ^ 6 when N — N b ~ Nf and ~ £/. In fact, for a repulsive boson-boson 
interaction, the boson distribution become broader than that for a free system just discussed, but this effect does not 
change the qualitative structure as shown in sec. V. On the other hand, the strong boson-fermion interaction may 
give qualitative changes of the density profiles such as the phase separation and the fermion collapse as shown in J7|]8| . 
We also discussed about them from different points of view in iii) in this section. 

ii) Scaled GP and TF equations and changes of the density distribution for the different N b j 

For qualitative estimations of the density distribution, we should use the scaled dimensionless variables || : 
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%=T> Pb,f = -Tf-nbJ, P>b,f — T u bJ, (9) 
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where N b j, £ b j are the total boson/fermion numbers and the harmonic-oscillator lengths defined in eqs. (|^) and (^). 
Using these scaled variables, the GP and TF equations in eqs. (|5|) and (||) become 

[-Vl + x 2 +gN b pb(x) + hN f pf{x)}ipb(x) = fi b ip b (x) (10) 

-j^[Qir 2 N f p f (x)}i +R m Rlx 2 + hN b pb{x) = p f , (11) 

where R m = mf/m b and R u — LOfjujb- The scaled coupling constants g and h in eqs. ( |l0| ) and ( |li"| ) are defined by 

~_ 2 9 Q a bb 1 _ 2/i a b f ( 1 \ 

where eq. (^) has been used to obtain the a^^-representations, and the last factor comes from mb/rribf = 1 + i?^, 1 . 

From eqs. ( |l0| ) and (|TT1), we should notice that the scaled equations include six independent parameters, (gNb, 
hNf, Nf, hNb, R m , Rw)\ for two systems where these parameters take the same values, their physical properties 
become similar under the scaling relations ([)]). 

Here, based on eqs. ( fio| ) and (|ll|), we discuss about changes of the density distribution by the scaled parameters. 
In this subsection, we assume R m = R^ = 1 and the large particle numbers N b ~ Nf ~ N 3> 1; for the effect of 
R w 7^ 1, the discussion will be given in sec. V. 

Let us consider the fermion TF equation (|ll|), first. For the chemical potential fif, we use the result in the system 
without boson-fermion interaction: p,f ~ 2(6Np ) 1 / 3 ~ TV 1 / 3 . It should be a good approximation when Nf ^> 1. Using 
it in eq. (|l0|), we obtain pf ~ NJ 1 ^ 2 at the central region (x ~ 0). On the other hand, the boson equation ( ^p| ) gives 

Ph ~ (gNb)~ 3 / 5 under the TF approximation valid in Nb S> 1. Using these estimates, we can obtain the ratio Rmt of 
the boson-boson/boson- fermion interaction effects in the GP equation (|l0|): 
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x=0 ~9 N b {gN b )-% gN b V ^/ 



(13) 



The small Ri nt (<SC 1) indicates the small contribution from the boson-fermion interaction for the boson part, and 
we obtain pb ~ Pb(/J = 0). However, when Ri n t > 1, the boson-fermion interaction is not negligible, and the boson 
density profile can be quite different from that in non-interacting case (h — 0). 
In 39 K- 40 K system discussed in sec. V, the scaled parameters are 

{gN b , hNf, Nf, hN b ) = (2.68 x lQ- 2 N b , 1.57 x 10~ 2 iV/, N f , 1.57 x 10~ 2 N b ). (14) 

In N ~ N b ~ A^, gA^ and ftJVf are in the same order, but Nf is a factor of 2 larger than WV/, so eq. ( |l3| ) gives 
•Rint ~ 10" 1 . It explains the rather small change of the boson density profile pb ~ p&(/i = 0) in 39 K- 40 K system, as 
shown in the numerical results in sec. V. 
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It is very interesting to apply the above estimation for the results by |Q|S|] and compare them with the results in 
this paper. One of the parameter set with which the boson density profiles changed drastically is 

{gN b , hN f , N f , hN b ) = (4.22 x 10 6 , 4.22 x I0 6 h/g, 10 6 , 4.22 x 10 6 h/g), (15) 
where h/g = h/g = ~ 5/4 has been taken and the qualitative changes of the boson density profile have been seen 



h/g > 3/4 |7jj8|]. It should be noticed that the parameters in (15) are in the same order when h/g ~ 1, and it means 
Rint ~ 1 in that case. For this reason, the boson density profile changed drastically from that of h — with these 
parameters. 

To check the relation between the boson profile change and the ratio Ri n t in eq. (|l^) a little more, we solved eqs. 
( fTo| ) and ( |lT|) numerically for a) Nf — 10 6 and b) Nf — 10 s with other parameters fixed as in (15), and, for h/g, two 



values were taken h/g = 0,0.6. The resultant boson densities are shown in fig. 1: n b (h ^ 0) has very different profile 
from n b (h = 0) in case a), but the very small change is found between them in case b). Those results shows the clear 
relation between Ri nt and the boson density profile. 

In summary, the parameter regions exist where the boson density profile changes very large or not, and these regions 
are discriminated by the value of Rint- 

We should comment that the results in this paper are complementary with those with (|l5|) in Q in the scale of the 
ratio R int . 

iii) Boson-fermion interaction effects for the fermion distribution function 

We next consider the effect of the boson-fermion interaction on the fermion distribution, especially the phase 
separation phenomena in the strongly repulsive boson-fermion interaction case, which has been proposed in Jt]||. Here 
we assume the TF approximation for bosons (neglecting the kinetic term in eq.(|])). As discussed in the literature 
[0 ^) , the effective potential for fermions receives an additional repulsive or attractive contribution from bosons within 
the range r < R b = y '2fi b / 'm^w 2 depending on the sign of the interaction. 

For a strongly repulsive boson-fermion interaction (h > 0), the fermion will be squeezed out from the center, so that 
the two kind of particles tend to make a separate phase |t],|| . Here we give analytical estimation for this interesting 
phenomena. As shown in , in the phase-separated system, the fermion are almost completely pushed away outside 
the boson distribution, so that the critical ratio (h/g) c for the phase separation can be estimated from the vanishing 
fermion density at the center n/(r = 0) = Q. Neglecting the boson-fermion interaction in the boson TF equation, we 
obtain 

N f = F(a) = ( fa(a-2)(a-l) 1 / 2 + a 3 arctan(a-l) 1/2 | (16) 

where a = h/g > 1. In fig. 2, the function F(a) in ( [l6| ) is plotted when Nf = 10 3 with the parameters used in ||]: 
(gN b , hNf, N f , hN b , R m , R u ) =(3.0 x 10 4 , 3.0 x 10a, 10 3 , 3.0 x 10a, 1, 1) @. From the crossing point of /(a) and 
Nf = 10 3 (dashed line) in fig. 2, the critical point can be read off: a = (h/g)\ c ~ 1.3, which is almost consistent with 
the result given in || (a ~ 1.18). 

For an attractive interaction (h < 0), fermions tend to concentrate and increase the overlap with bosons. Thus we 
may expect a coherence of the two kinds of particles to occur in this case. For large attraction, another interesting 
phenomena, fermion collapse in mixture has been proposed M. 



IV. NUMERICAL PROCEDURE 



In the numerical calculation, we solved the set of equations (|^) and (^|) by reducing them into the following equation 
for the order parameter: 



d 
~dt 



J it? 

9—r- 



24tt 2 



m h io\ 



3/2 



f(t) = o, 



(17) 



where t = x 2 = (r/£b) 2 and the scaled parameters x, g, h and fibj have been defined in (j^). The function / is related 
to the order parameter by $(r)/r = e~ x f(x 2 ). The boundary conditions for f(t) are given by the value at t = 
and by the asymptotic condition f(t — > oo) ~ ^CAb— 3)/4 ^ e jj rg ^ gj ve m ^j a i values for /(0) and fif and solved the 
equation ( p"7| ) numerically by means of the relaxation method. Finally, the fermion density nf(r) is obtained from 
eq.(H). The details of the calculation will be given elsewhere p5[. 



4 



We consider the potassium atoms where there exist two bosonic ( 39 K, 41 K) and one fermionic ( 40 K) isotopes. 
The precise values of the interatomic interaction are not known, but the recent estimate from molecular scattering 



L2 13 1 suggest a repulsive interaction between 39 K and 40 K, and an attractive one between 41 K and 40 K, although 
the values still have rather large ambiguities. We take up the following values for the s-wave scattering lengths fl3| : 
a 66 ( 39 K- 39 K) = 4.29nm, a fc/ ( 39 K- 40 K) = 2.51 nm, a 66 ( 41 K- 41 K) = 15.13 urn, a h/ ( 41 K- 40 K) = -8.57 nm. We take 
the atomic masses to be the same for all the isotopes: m b = mj = 0.649 x 10~ 25 kg. The angular frequency w b of 
the bosonic external potential is fixed at 100 Hz, while w/ is allowed to vary. In this case, it is corresponding to 
g = 2.68 x 10~ 2 , h = 1.57 x 10~ 2 for 39 K- 40 K system and g = 9.44 x 10~ 2 , h = -5.34 x 10~ 2 for 41 K- 40 K system. 

V. RESULTS 

First, we consider the 39 K- 40 K system where the boson-fermion interaction is repulsive. In this case, because of 
hN b /gN b < 1, phase separation do not occur. Fig. 3 shows the density-distribution functions of bosons (a) and 
fermions (b) for A^ = 10000, Nf = 1000 and u>f = oj b = 100 Hz. The dashed line shows the result for a b f = 
(h = 0) as compared with the solid line for a b f ^ (h ^ 0). As seen from the figure, the fermions have a much 
broader distribution than bosons even for a much less number of particles. That has been discussed in sec. Ill-i). The 
boson-fermion interaction is seen to squeeze out the fermions leading to a fermion-density depletion at the center, 
while its effect on the boson distribution is negligible. This also can be seen in Table I where some observables are 
listed for the two cases (N b , AT/) = (1000,1000) and (10000,1000) with/without the boson-fermion interaction. The 
robust structure of the boson distribution just mentioned is reflected in the almost constant values of the rms radius 
and other observables for the bosonic part of the mixed condensate against the switching on/off of the boson-fermion 
interaction. To understand this robustness, we apply the estimation in sec. Ill-ii): the (1000,1000) system just 
corresponds to the parameters of eq. (|l4|), so the Ri nt in ( |l3| ) becomes Rint ~ 10~ 4 . For (10000,1000) case, the scaled 
parameters are (gN b , hNf, Nf, hN b ) =(2.68 x 10 2 , 1.57 x 10, 10 3 , 1.57 x 10 2 ), and the ratio hN s /gN b in eq. (]l|) gives 
additional factor 10 _1 , so we obtain Ri nt ~ 10~ 2 . It shows that (10000,1000) system is more robust than (1000,1000) 
case. 

It is supported in ref. || that the TF approximation yields qualitative correct result. Using the method in |^6| , we 
can also verify that the correction for the TF approximation is small for 1000 fermion case. 

We next turn to the 41 K- 40 K system where the boson-fermion interaction is attractive. Because of the large value 
of the strength, one may expect a sizable effect of the boson-fermion interaction for this system. Moreover, the 
strong boson-boson interaction makes a stronger overlap of bosons and fermions and thus the attractive interaction 
are more efficient. So the system might be expected to be unstable against collapse, if the effect of the boson-fermion 
interaction is too strong. Actually, (N b , A r /)=(1000,1000) and (10000,1000), we obtained a distribution as given in 
fig. 4 and Table II. The effect of the attraction on the fermions is larger than the 39 K- 40 K system, although it is 
not strong enough to cause an appreciable change in the boson distribution. The stability condition for the fermion 
collapse can be derived from the scaling law in refs. 0: 

— < c = 2 x 6.9AT- 1 / 6 . (18) 
9 



In the present case, we obtain h 2 /g = 3.02 x 10~ 2 and c ~ 4.4, so that the stability condition ( |l8| ) is well satisfied, 
and we have found the parameters of 41 K- 40 K system is in (quasi) stable region. 

Is the bosonic distribution always robust against a mixture of fermions? Aside from the possibility of changing the 
interaction strength via Feshbach resonances fTij] , we have a possibility of enhancing the boson-fermion interaction 
effect using the different confining potential for bosons and fermions separately (i.e. LOf/oj b 7^ 1). Fig. 5 shows the 
boson/fermion distribution for (a) 39 K- 40 K (R u = 1,5,10) and (b) 41 K- 40 K (R u = 1,4,7): {N b ,N f ) = (1000,1000), 
uj b = 100 Hz. For a large value of the ratio R u — u>f/u> b , the fermions are trapped in a strong confining potential and 
thus have a large overlap with bosons; the boson distribution (solid lines) is sensitive to the value of R^ even though 
the bosonic parameters are fixed. This may be traced back to a rather large fraction of the boson-fermion interaction 
energy in the total boson energies: for the 41 K- 40 K system at R u = 7, for instance, the ratio of the boson-fermion 
interaction energy to the total boson energy becomes Ebf /E to t( 41 K) — 0.36 |Q, while it is less than 1 per cent to 
the total fermion energy. The change of the boson-density distribution can be estimated from the ratio Rint in Jl3]). 
When R u ^ 1, it becomes 

hN f f~gN b \i 3/2 
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so that the boson-density distribution become more sensitive for large values of R^, . 

The effect of the boson- fermion interaction on fermions may show up in the opposite extreme, i.e., for vanishing 
R u . In this case, the fermions can keep themselves from escape only through the attractive boson-fermion interaction. 
This possibility can be estimated from the chemical potential: when fermions are trapped by bosons, the [if should 
take the negative value. Taking the TF approximation for the bosonic part (gNf, 1, gNf, 3> hNf) and setting 
ujf = for the fermionic part, we obtain the total fermion number represented by [if. 



N f = 



1 



48\h/g\ 



J 1 M/ 



2/5' 



(20) 



When [if = in (|20|), we obtain the maximum number for the boson-trapped fermions: 



N max 



1 

48 



3/2/,, x 6/5 



2.8 x 10- 4 7V b 6/5 



(21) 



where we used the parameters of 41 K- 40 K system to obtain the last term. From eq. (|2l]), estimations can be obtained 
for the boson-trapped fermion number: for Nf, = 10000, only 18 fermions are bound, while for TV;, = 10 6 more than 
4000 fermions will be trapped within the bose condensate. This qualitative estimate is verified also in the numerical 
calculation for Nf, — 10000, although one should actually take the shell effect into account for such a small number of 
fermions [fill. 



VI. SUMMARY 



In the present paper, we studied static properties of a mixed system of bosons ( 39 K or 41 K) and fermions ( 40 K) in 
the trapping potential at temperature zero. We solved a coupled GP and TF equations for a different combination 
of boson and fermion particle numbers and also for a repulsive or an attractive interaction of bosons and fermions 
depending on the combination of potassium isotopes. 

The analytical estimations have been done generally for the rms radii of the boson/fermion distribution, the 
parameter-dependence of their profile and the phase separation, and we have found the feasible conditions for them. 
These analytical estimations have been applied to the numerical results for the potassium system and also compared 
with the discussion done in [0J^]: consequently, the consistency was checked for these estimations. 

For similar particle numbers and external potential parameters, the fermions have a much more extended distribu- 
tion and larger energies than bosons because of the exclusion principle. We studied the effect of the boson-fermion 
interaction on the density distribution of both kinds of particles. For a repulsive interaction, the fermions tend to be 
pushed out, but the effect is not very large for the parameters of 41 K- 40 K system adopted in the present paper. In the 
case of the attractive interaction, a coherence of 41 K and 40 K occurs and is strongly enhanced if the trapping potential 
is adjusted so as to make them overlap. It was shown also that the fermions could be kept trapped without external 
potential due only to an attractive interaction with trapped bose particles. One may as well consider the opposite 
case where the bosons are trapped only via an attractive interaction with surrounding trapped fermi particles. This 
will be discussed in a future publication 
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TABLE I. Results for the 39 K- 40 K system with (a b f = 2.51 nm) and without (a b f = 0) the boson-fermion interaction 
for (Nt,Nf) — (1000, 1000) and (10000, 1000). Contributions to the equilibrium energy (measured in the unit of hu> b /2) per 
particle, the central density n(0) (x = r/£ b ) and the rms radii are given for each of the constituent isotopes. Other parameters 
are fixed as abb = 4.29 nm, LOf — co b = 100 Hz. 



(N b ,N f ; a b f) 


atom 


fib or /j/ 


Etot/N 


Ekin/N 


E ho /N 


E BB /N 


Ebf/N 


n(x = 0) 




(10 3 ,10 3 ;0.00) 


39 K 


4.19 


3.66 


1.16 


1.96 


0.53 


0.00 


98.7 


1.40 




40 R 


36.3 


27.3 


13.6 


13.6 




0.00 


3.70 


3.69 


(10 3 ,10 3 ;2.51) 


39 R 


4.24 


3.71 


1.16 


1.96 


0.53 


0.05 


98.7 


1.40 




40 K 


36.4 


27.3 


13.6 


13.7 




0.05 


3.47 


3.70 


(10 4 ,10 3 ;0.00) 


3y K 


8.21 


6.26 


0.70 


3.62 


1.95 


0.00 


291 


1.90 




40 R 


36.3 


27.3 


13.6 


13.6 




0.00 


3.70 


3.69 


(10 4 ,10 3 ;2.51) 


39 K 


8.25 


6.31 


0.70 


3.62 


1.94 


0.05 


291 


1.90 




40 K 


36.6 


27.7 


13.3 


14.0 




0.45 


3.07 


3.74 


TABLE II. Results for 


the 41 K- 40 K 


system. 


See Table I 


for the notation. The parameters 


are: abb = 


15.13 nm, 


u)f = (Jb = 100 Hz. 




















(N b ,N f ; a bf ) 


atom 


fib or fi f 


Etot/N 


Ekin/N 


E ho /N 


Ebb/N 


Ebf/N 


n(x = 0) 




(10 3 ,10 3 ;0.00) 


41 K 


5.83 


4.68 


0.90 


2.62 


1.15 


0.00 


53.5 


1.62 




40 R 


36.3 


27.3 


13.6 


13.6 




0.00 


3.70 


3.69 


(10 3 ,10 3 ;-8.57) 


41 K 


5.65 


4.49 


0.91 


2.62 


1.16 


-0.19 


53.9 


1.62 




40 R 


36.2 


27.1 


13.8 


13.5 




-0.19 


4.13 


3.68 


(10 4 ,10 3 ;0.00) 


41 K 


13.0 


9.58 


0.50 


5.65 


3.43 


0.00 


135 


2.38 




40 K 


36.3 


27.3 


13.6 


13.6 




0.00 


3.70 


3.69 


(10 4 ,10 3 ;-8.57) 


41 K 


12.9 


9.40 


0.50 


5.62 


3.46 


-0.18 


136 


2.37 




40 R 


35.4 


25.6 


14.7 


12.7 




-1.78 


4.71 


3.56 



FIG. 1. Changes of the boson density distribution for different Nf values: a) Nf = 10 6 and b) Nf = 10 8 with (gNb, hNj, 
hN b , R m , R„) = (4.22 x 10 6 , 4.22 x 10 6 |, 10 6 ,4.22 x 10 6 |, 1, 1). Solid lines are for h/g = 0.6 and dotted lines are for h/g = 
(no boson-fermion interaction). 



FIG. 2. Critical line of component separation for the system with (gN b , hNf, Nf, hN b , R m , Ruj) =(3.0 x 10 4 , 3.0 x lQh/g, 
10 3 , 3.0 x IQh/g, 1, 1). 



FIG. 3. Density distribution for the ( 39 K, 40 K) = (10000,1000) system: (a) for 39 K and (b) for 40 K (b) are plotted against 
x = r/^b- The solid line shows the full calculation with the (repulsive) boson-fermion interaction (a b f = 2.51 nm) and the 
dashed line without. Both lines are indistinguishable for the bosonic part (a). 

FIG. 4. Density distribution for the ( 41 K, 40 K): (a) for the 41 K and (b) for 40 K. The boson-fermion interaction is attractive 
in this case (a b f = — 8.57nm). 

FIG. 5. Dependence of the density distribution on the oscillator- frequency ratio R^ — LOf/cjb- a) 39 K- 40 K system with 
R^ = 1,5, 10 and b) 41 K- 40 K system with i?„ = 1,4,7, where (N b ,N f ) = (1000, 1000) and cu b = 100 Hz. The solid line shows 
the boson density n b (x) and the dashed line shows the fermion density rif(x). 
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